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\q ' Abstract 

, We discuss the initial-boundary value problem of General Relativity. 

Previous considerations for a toy model problem in electrodynamics mo- 
tivate the introduction of a variational principle for the lapse with several 
attractive properties. In particular, it is argued that the resulting ellip- 
tic gauge condition for the lapse together with a suitable condition for 
the shift and constraint-preserving boundary conditions controlling the 
Weyl scalar are expected to yield a well posed initial-boundary value 
problem for metric formulations of Einstein's field equations which are 
commonly used in numerical relativity. 
■ To present a simple and explicit example we consider the 3 + 1 de- 

composition introduced by York of the field equations on a cubic domain 
with two periodic directions and prove in the weak field limit that our 
gauge condition for the lapse and our boundary conditions lead to a well 
posed problem. The method discussed here is quite general and should 
also yield well posed problems for different ways of writing the evolution 
equations, including first order symmetric hyperbolic or mixed first-order 
second-order formulations. Well posed initial-boundary value formula- 
tions for the linearization about arbitrary stationary configurations will 
be presented elsewhere. 
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1 Introduction 

A common approach to numerical relativity is to find solutions of Einstein's 
field equations on a cylinder of the form / x X, where / = [0, T] is a time 
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interval and £ is a three-dimensional compact manifold with smooth boundary 
9S. The initial-boundary value problem (IBVP) consists in finding solutions 
with specified (initial) data on the "bottom" of this cylinder (i.e. on {0} x S) 
and boundary data on the set T = I x 9S. The exact form of the boundary 
data depends on the nature of T; in this article we consider the case where T is 
timelike. This excludes the case of inner black hole excision boundary conditions 
where T is spacelike or null, and the case where spacetime is compactified in 
such a way that T is a null surface. 

For timelikc T, boundary conditions must be specified with care. First, they 
should be stable in the mathematical sense of well poscdncss of the resulting 
IBVP which means that one can guarantee local in time existence of solutions 
which lie in an appropriate normed function space and which depend uniquely 
and continuously on the initial and boundary data. Second, the boundary con- 
ditions should be constraint-preserving which means that they must guarantee 
that a solution to the IBVP with constraint-satisfying initial data automati- 
cally satisfies the constraints everywhere on / x S. Third, the boundary data 
should enable control over the gravitational flux through the boundary, at least 
in some approximate sense. This is a difficult question since in General Rela- 
tivity there are no meaningful local expressions for quantities representing the 
energy density of the gravitational field. This is related to the fact that locally, 
the gravitational field can always be transformed away. Finally, it would be 
desirable that boundary conditions could be specified in such a way to control 
part of the geometry of T. 

A well posed IBVP of Einstein's vacuum field equations addressing all the 
above points has been presented in pQ . This work is based on the use of tetrad 
fields in a gauge adapted to the boundary T and casting the evolution equa- 
tions into a first order symmetric hyperbolic system for the tetrad fields, the 
connection coefficients and the components of the Weyl tensor. Furthermore, 
the evolution equations in PP imply that the constraint propagation system, de- 
scribing the evolution of the constraint fields, has the form of a hyperbolic system 
which on T contains only derivatives tangential to T. As a consequence, the 
satisfaction of the constraint-preserving property is automatic. The only fields 
with non-trivial characteristic speeds at the boundary are the Weyl scalars *o 
and ^4 that can be interpreted as describing in- and outgoing gravitational 
radiation, at least in some approximate sense. (Notice that our definition of 

and ^4 differs from the one in pQ.) By specifying data to or a suitable 
combination of Wo and ^4, the work in Ref. pQ obtains a well posed IBVP for 
Einstein's vacuum field equations by applying standard theorems on symmetric 
hyperbolic equations with maximally dissipative boundary conditions |2J 03 0] - 

On the other hand, most numerical codes are based on formulations of 
Einstein's equations which use metric variables and not tetrad fields. For 
such formulations, the problem of obtaining a well posed IBVP is open. One 
of the major obstacles stems from the fact that compatibility with the con- 
straints yields boundary conditions which have the form of differential equations 
on the boundary for which no standard mathematical theorems are known. 
Despite several partial results El El El HO EH]- the derivation of necessary 
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conditions for well posedness [111 1121 1131 |§1 an d numerical experiments 
HH El El 113 113 UHl C2l EDI El , a full understanding of the IB VP for metric 
formulations has not yet been achieved. 

In this article, we reconsider the IBVP for Einstein's field equations by com- 
bining the ideas of Ref. pp with new ideas developed in [5T] for a model problem 
in electrodynamics. These new ideas, which are reviewed in Sect. |3 are based 
on separate estimates for the constraint fields, the curvature fields and the im- 
position of an appropriate gauge condition which allows an estimate for the 
remaining, gauge-dependent fields. Boundary conditions on the constraint and 
curvature fields are an essential part in obtaining these estimates. In the electro- 
magnetic case, the curvature fields are the electric and magnetic fields which are 
gauge-invariant. In General Relativity, the curvature field is described by the 
Riemann tensor which can be decomposed into the Ricci tensor and the Weyl 
tensor. Motivated by the model problem in j5J we consider a gauge condition 
for the lapse that is obtained by minimizing a simple functional representing the 
norm squared of the time derivative of the extrinsic curvature. The minimiza- 
tion principle, which yields a fourth order elliptic boundary value problem for 
the lapse, possesses several useful properties. One of them is that the functional 
is positive semi-definite and zero if and only if the components of the extrinsic 
curvature are time- independent. Therefore, the gauge condition should be well 
adapted to stationary configurations or small deviations thereof. In Sect. 12.61 
we summarize the boundary conditions, the gauge conditions and the evolution 
system considered in this article and write down the resulting IBVP. In Sect. El 
we examine the variational principle for the lapse in the weak field limit, estab- 
lish its well posedness in Theorem ^ and give a physical interpretation for it as 
a gauge-fixing procedure for the extrinsic curvature. In Lemman we also derive 
some key estimates for the extrinsic curvature which are valid in this gauge. 

The main result of this article is derived in Sect. 0] There, we consider 
the linearized Einstein vacuum field equations in the 3 + 1 decomposition as 
described in Ref. [22] and couple them to the fourth order elliptic gauge condi- 
tion for the lapse. We prove in Theorem [3| that this elliptic- hyperbolic system, 
together with constraint-preserving boundary conditions controlling a suitable 
combination of the Weyl scalars an d ^4 leads to a well posed IBVP. The 
proof is based on casting the evolution system into an abstract Cauchy problem 
on an appropriate Hilbcrt space and verifying that the evolution vector field 
generates a strongly continuous semigroup. Possible generalizations to other 
formulations of the linearized Einstein equations and future work are discussed 
in the conclusions, Sect. |3J Appendix ^ summarizes standard results from 
Fredholm theory used in the discussion of the gauge condition for the lapse, and 
Appendix IB1 contains technical proofs of some of the statements made in Sect. 

m 
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2 Main ideas 



In this section we sketch some ideas for obtaining a priori estimates in the IB VP 
of Einstein's field equations. We start in Sect. 12.11 bv reviewing the electromag- 
netic toy model problem analyzed in |21| and showing that the gauge condition 
used in that work can be obtained by variation of a simple functional. In Sect. 
12. 21 we decompose Einstein's field equations into evolution and constraint equa- 
tions. Next, in Sect. 12.31 the constraint propagation system describing the 
evolution of the constraint fields is derived and a boundary condition which 
yields estimates for these fields is constructed. In Sect. 12.41 we discuss how 
estimates for the curvature fields can be obtained. The curvature tensor can be 
decomposed into the Ricci tensor and the Weyl tensor. The Ricci tensor can be 
estimated through Einstein's equations provided suitable estimates on the mat- 
ter fields are available. Therefore, it remains to estimate the Weyl tensor. In 
order to do so, we analyze its propagation and explain the idea for constructing 
boundary conditions and deriving estimates for the Weyl tensor following the 
lines of Ref. pQ. Next, in Sect. 12.51 the first and second fundamental forms 
are estimated. For this, an appropriate gauge condition for the lapse is needed. 
Such a condition will be derived by introducing a functional in analogy with the 
electromagnetic case. Finally, the resulting IBVP for Einstein's vacuum field 
equations is summarized in Sect. 12.61 The well posedness of this problem in the 
weak field limit is proven in the subsequent sections. 

In the following, we assume that £ is an open bounded subset of R 3 with 
C°° boundary <9E. The spacetime metric on [0, T] x £ has signature (— , +, +, +) 
and is denoted by g and the unique metric-compatible, torsion-free connection 
by V. The induced three-metric and connection on X are denoted by h and D, 
respectively. The indices a, 6, c, d, e, / denote spacetime indices running from 
to 3 while i, j, k, I denote spatial indices running from 1 to 3. We shall use the 
notation = for equalities holding on 9£ or on [0, T] x 9£. The unit outward 
normal one-form to <9X is denoted by iVj and = 8^ — NiN? is the projector 
onto the tangent space of dH. 

2.1 The electromagnetic case 

In |21) Maxwell's equations are written as a first order system for the magnetic 
potential Aj, the electric field Ej and the derivatives of the magnetic potential, 
Wij = DiAj. In the source-free case, this system is 



subject to the constraints C = D l Ei = and Cy = Wij — D^Aj = 0. Here, <f> 
denotes the electrostatic potential, hij the components of the Euclidean metric 
and a and r are nonvanishing constants having the same sign. 



dtA. 
d t E L 

d t Wi, 



Ej + Dj<j>, 

D l W t] ~ (1 + tr)DW J( + ah M DjW u , 
DiEj + T - h lJ D k E k + DiDrf, 



(1) 
(2) 

(3) 
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A well posed IB VP for this system is obtained in the following way: First, the 
constraint propagation system describing the time evolution of the constraint 
fields C and dj is analyzed. It is given by 

StCn = \hijC, (4) 

d t C = -aD'Q, (5) 
d t d = -rDiC, (6) 

where Cj = D^Cij — h kl DiCki- By imposing, for instance, the boundary con- 
dition N l Ci = one obtains an L 2 estimate for the constraint fields. Next, 
using this, the physical energy of the system and the radiative-type boundary 
condition 

-r^Ei-N^Wa-W^^gi, (7) 

where gi is a given vector field on the boundary satisfying N l gi = 0, one obtains 
L 2 bounds for the electric and magnetic fields, Ei and Wuj] , respectively. The 
remaining difficulty is to estimate the symmetric, gauge- dependent, part of Wij . 
It turns out that the simplest gauge choice, the temporal gauge <j> = 0, does not 
lead to such an estimate. This can be seen by means of the following family of 
electrostatic solutions [21] 

A j = tDjf, 

Ej = Djf, (8) 
Wij = tDiDjf, 

where / is a smooth, time- independent, harmonic function. This family solves 
the evolution and constraint equations and the radiative-type boundary condi- 
tion with boundary data gi = j^Djf. A problem is the time-dependence of 
the magnetic potential which implies that for t > the solution depends on 
second derivatives of / while the initial and boundary data depends only on 
first derivatives of /. This violates the expected energy estimate in L 2 . 

For this reason, in |21j. a different gauge condition was introduced which 
can be obtained by minimizing the functional 

over the space H 1 ^) consisting of square- integrable functions <f> on E which 
have square-integrable first order spatial derivatives. Therefore, this gauge con- 
dition minimizes the L 2 norm of the time-deformation of the magnetic potential 
over all possible electrostatic potentials <j) for which this norm is defined. Here 
and in the following, \fh and denote, respectively, the determinant and 
the inverse of the three-metric. In the above example, where Ej = Djf, the 
minimum of I em is zero and yields (f> = —f which in turn implies that Aj is 
time- independent. Therefore, this gauge condition precludes the above coun- 
terexample. Furthermore, it was shown in |21| to yield a well posed IBVP. A 
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straightforward generalization of Iem is the family of functionals 



Iem, 



M = l I [n 2 tf +h l J(d t A l )(d t A 3 )]Vhd 3 x 



parametrized by a nonnegative constant /i. The extra term, /i 2 </> 2 , in the inte- 
grand plays the role of a penalty term which prevents \(/)\ from becoming "too 
large" and might be useful for numerical simulations. Variation of Iem,^ with 
respect to <p yields the elliptic boundary value problem 



Notice that for p, > the solution of this boundary value problem is unique 
while for /i = the solution <j) is only unique up to an additive constant which 
might cause problems in an elliptic solver. The important point here is the 
boundary condition (|10|l which allows to set the normal component of the mag- 
netic potential to zero at the boundary: Aj = 0. This boundary condition 
in turn allows to estimate the L 2 norm of A and its first order spatial deriva- 
tives in terms of the divergence and the curl of A. The divergence of A can be 
estimated by taking the divergence of Eq. JQ) and using Eq. © . The curl of A 
is known since it corresponds to the magnetic field Wuj] for which an estimate 
has been obtained. 

Using the methods of Ref. one can show that the IBVP described by 
the evolution equations Ill2l3fl for ctt > 0, the boundary conditions @ and 
N l d = and the gauge condition (|9I1()|I is well posed. The interesting point is 
that one still obtains a well posed problem in the limiting case a = t = where 
the evolution equations are weakly hyperbolic 1231 . In this case, the constraints 
propagate tangentially to the boundary, and the boundary condition N l d = 
has to be dropped. 

2.2 Einstein's field equations 

Next, we discuss how the ideas above can be applied to the IBVP of General 
Relativity. In the remaining part of this section, we consider the full nonlin- 
ear Einstein equations with matter fields and draw particular attention to the 
propagation of the constraint fields, the propagation of the Weyl tensor and the 
gauge conditions. An IBVP for the nonlinear case is mentioned at the end of 
this section. In the next two sections, we prove that this problem is well posed 
in the weak field regime. 

We start with Einstein's field equations in the 3 + 1 formulation of Ref. 
|22j . This is done for simplicity; as discussed in the conclusions, other 3 + 
1 formulations might serve as a starting point. Assume spacetime (M,g) is 
globally hyperbolic, i.e. there exists a globally defined time function t ; M —> I 
such that M is foliated by spacelike hypersurfaces E T = {p £ M : t(p) = r} . Let 
n a = —aSJat be the future-pointing unit normal to these slices (we choose the 



+ D l Di)<t> + D i Ei = 
N\D t (f, + E l ) = Q 



on S, 
on <9E. 



(9) 
(10) 



G 



time orientation such that the lapse function, a, is strictly positive). The three- 
metric h a b and extrinsic curvature k a b (first and second fundamental forms) are 
defined by 

h a b = 9ab + n a n b , (11) 
kab = -V a n fc - n a a b , (12) 

respectively, where a& = n a V ' a n b = £>b(loga) is the acceleration field. The 3+1 
form of Einstein's equations yield the evolution equations 

£ n h a b = -2k ab , (13) 

C n kab = Rab D a Dba + kkab - 2k a d kbd - ns a b , (14) 

a 



and the constraints 



H = ^(h ab R ab + k 2 ~k ab k ab )- K p = 0, (15) 

M a = D a k - D b k ab + KJ a = 0, (16) 



where here and in the following, D a and R a b denote, respectively, the covariant 
derivative and the Ricci tensor compatible with the three-metric h a b, k = h ab k a b 
is the trace of the extrinsic curvature, k — 8irG where G is the gravitational 
constant, and where in terms of the stress-energy tensor T a b the source terms 
p, j a and Sab are given by 

a brri 

p = n n Tab , 
ja = -h a c n d T cd , 

Sab = (ha C hb d -h ab 9 Cd )T cd . 

For consistency with the Bianchi identities, the stress-energy tensor must be 
divergence-free which translates into the conditions 

C n p + ^D a (a 2 Ja ) - 2kp - (k ab - kh ab ) Sab - 0, (17) 

Cnja H ^D a (a 2 p) - kj a + — D b (as ab - ah ab h cd s cd ) = 0. (18) 

or a 

A well-known property of the evolution system (|13I14J) is that the Ricci opera- 
tor, h a b Rab, is not elliptic. In particular, this implies that simple algebraic 
gauge choices yield a weakly hyperbolic system for which there are no general 
well posedness results and standard discretizations lead to unstable numerical 
schemes |21]. A possible remedy to this problem is to adopt differ- 

ent gauge conditions, like the full harmonic gauge |27| . in which the evolution 
equations become manifestly hyperbolic. 
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2.3 Estimates for the constraint fields 



Independent on what gauge is chosen, the evolution system (|13I14|) has the 
property that it induces "nice" evolution equations for the constraint fields: As 
a consequence of the (twice contracted) Bianchi identity, the evolution equa- 
tions (|13I14|> and the consistency conditions (|17I18|I one finds that H and M a 
propagate according to the simple symmetric hyperbolic system 



which closely resembles a wave equation written in first order form. By speci- 
fying suitable boundary conditions for this system, estimates for the constraint 
fields H and M a can be obtained. In the following, we shall assume that the 
timelike boundary T is such that n a v a > everywhere on T, where v a denotes 
the outward unit normal 1 to T . As can be seen by a simple energy argument 
(cf. Lemma |21 below), the boundary condition 



guarantees that H = and M a = everywhere on [0, T] x £ provided H = 
and M a = on the initial surface {0} x S. Furthermore, the boundary condition 
(|21|) yields an L? estimate for the constraint fields H and M a . As we will see 
shortly, it turns out that we need more regularity and also need an L 2 estimate 
for the first spatial derivatives of H and M a . It will be shown in the linearized 
case that this can be achieved using the boundary condition (|21|) . 

Before we proceed, we make two comments. First, we stress that the evo- 
lution equations derived by Arnowitt-Deser-Misner 2S : differ from the one dis- 
cussed here by the addition of the Hamiltonian constraint to the evolution equa- 
tions for k a b- As a consequence, in the former case, the resulting constraint 
propagation system is only weakly hyperbolic which makes it much more dif- 
ficult to control the constraint fields. (See Ref. [2SJ for the details on this 
important difference.) Second, we would like to mention that if n a v a = and 
in the absence of matter fields, the boundary condition H21J1 corresponds to one 
of the four Einstein boundary conditions v a G a b = 0, where G a b denotes the 
Einstein tensor, proposed in Ref. jHOj. (In vacuum Eq. I|21|) is equivalent to 
v a n b G a b = 0.) However, the imposition of all four Einstein boundary conditions 
can lead to difficulties. For example, imposing the two conditions v a n h G a b = 
and v a v h G a b = overdetermines the constraint propagation system since by 
virtue of the evolution equations the second condition is equivalent to setting 
H = 0. Therefore, if the initial data violates the constraints (as is the case in 
most numerical simulations), a solution to the IB VP might not even exist. 

1 Different boundary conditions have to be specified for the case where n a v a < everywhere 
on T . We do not analyze this case here. 




(19) 



(20) 



v a M a = 



(21) 
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2.4 Estimates for the Weyl curvature 

The next step is to obtain estimates for the Weyl curvature tensor C abcd . In 
terms of the electric and magnetic part of the Weyl tensor, defined by 

Ebd = C abcd n a n c , (22) 
B bd = \n a C aheS ^ d , (23) 

where e bcd = n a e abcd denotes the natural volume element on (S t ,/i a b), the 
Bianchi identities yield 

C n E ab = -e cd{a (D c + 2a c )B d b) 

- 3k {a d E b)d - 2k d {a E b)d + 2kE ab + h ab k cd E cd + P ab , (24) 
C- n B ab = +e c d(a(D c + 2a c )E d b - ) 

- 3k (a d B b)d - 2k d [a B b)d + 2kB ab + h ab k cd B cd + Q ab , (25) 



where 



Pab = ^ [ D (* M b) + 2a(aM 6 ) + 2k ab H 



-k {C n s ab + k( a d s b ) d + k ab h cd s cd )] (26) 

1 

2 



Qab = -^[zcd{ak c b)M d + ne cd(a D c s d b) ], (27) 



and the constraints on E and B 



D b E ab + 2k cd e b da B cb = P a , (28) 
D b B ab - 2k cd e b da E cb = Q a , (29) 



where 



Pa ^ \D a H + \ 



k a b M b - kM a - k (D b s ab - h cd D a s cd 



(30) 



Q a = -^e a cd [D c M d + Kk c b s db \. (31) 

In the above expressions for P a , Q a , P ab and Q ab we have used the constraint 
propagation equations (|19I2Q(I in order to eliminate the Lie derivatives of the 
constraint fields H and M a with respect to n a . 

Eqs. I|24I25|I constitute a symmetric hyperbolic system for the fields E and 
B which bears a striking resemblance to Maxwell's equations. Notice that the 
source terms P ab and Q ab depend on the constraint fields H and M a and their 
first spatial derivatives as well as on first derivatives of the matter terms s ab . 
This is the reason why more regularity is needed for the constraint fields (and the 
matter fields). A well posed IBVP for the subsystem H24I25|) can be constructed 
using the method in Ref. £Q where suitable combinations of the constraint 
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equations are added to the evolution equations (|24I25[) (see also Sect. 14. 2|) . 
Boundary conditions for the resulting linear evolution system for E and B which 
lead to a well posed problem in L 2 can then be specified in the following form 
HJ: Set N a = h ai Ni and choose for each t G [0,T] two unit vector fields e% 
and 63 on {t} x <9£ which are mutually orthogonal, such that N a , e%, e% form 
an orthonormal set. Define m a = 2~ 1 / 2 (e2 + «e|) and introduce the complex 
Newman-Penrose scalars 

*o = [E ab - e cda N c B d b ] m a m\ (32) 
1-4 = [Eab + e cda N c B d b ] fh a m b . (33) 

Next, let c be a smooth complex- valued function on [0,T] x (9E with magnitude 
smaller or equal than one, and let q be a smooth complex-valued function on 
[0, T] x <9£. Boundary conditions can then be specified in the form [Tlll3j 

*o = c* 4 + q, (34) 

where an overbar denotes complex conjugation. The appearance of the com- 
plex conjugation in Eq. (|34|l can be understood as follows: Under a rotation 
777 1 ► e l<fi m of m, * >-> e 2 ^* while * 4 i-> e" 2 ^* 4 . Therefore, if q h-> e 2itfi q, 
the boundary condition lf3"4"j) is invariant with respect to such rotations. Asymp- 
totically, and ^4 represent, respectively, the amount of in- and outgoing ra- 
diation. Therefore, the function c can be interpreted as a reflection coefficient. 



2.5 Estimates for the first and second fundamental forms 



So far, we have obtained estimates for the constraint fields and the curvature 
fields E a b and B ab . It remains to estimate the three- metric h ao and the ex- 
trinsic curvature k a b- For this, we first notice that the electric and magnetic 
components of the Weyl tensor can be expressed as 



E a b 
B a b 



^cd(a 



^rikab 

D c k d 



Rab 



1 



-D a D b a + kk at 



TF 



b) 



(35) 
(36) 



where [■■■} TF denotes the trace- free part with respect to h a b- This allows us to 
rewrite the evolution system (|13I14|) as 



£"n,fcab 



2,k n 



E„ 



1 



-D a D b a - k a c k c 



-hnhH 



^rSab + ^Kb (4p - h cd s cd ) 
2 



(37) 



(38) 



Since estimates for E ab and H have already been obtained, and provided suitable 
estimates for the matter terms p and s ab are available, one should therefore 
be able to integrate this system and obtain estimates for h ab and k a b- The 
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question here is what precise estimates we need. Since one would like to be able 
to estimate the Christoffel symbols, we need to control at least the L 2 norm 
of the three-metric and its first order spatial derivatives. Such an estimate can 
be obtained for fixed shift by integrating Eq. <|3*7|) provided a similar estimate 
for k a b is available. As in the electromagnetic case, one does not expect such 
an estimate to exist for any gauge condition for the lapse. Comparing Eq. 
with Eq. (|38|l and making the analogy A4 <-> k a b, £Jj <-> E a b, 4> <^r a between 
the electromagnetic case and General Relativity, the functional Iem suggests 
the consideration of the following functional over the space H 2 (Y?) consisting of 
square-integrable functions a on £ with square-integrable first and second order 
spatial derivatives, 

i GR [a] = y^h bd (^y^jVhd 3 x 

= \ f h ac h bd (-D a D b a - -C p k ab - F ab 
2 Jt. \a Gi 



s 

-D c D d a - -Cpk cd - F cd ) Vhd 3 x, 
a a 

where F ab = R ab + kk ab — 2k a d k b d — Ks ab and f3 a is the shift vector field. The 
reason for introducing the factor X/a in front of dtk ab is that without this factor, 
the minimum of Iqr is always zero and is attained for a — when Cfjk ab van- 
ishes 2 . Notice that minima need not be unique: For example, if Lpk ab = 0, the 
functional is invariant with respect to the rescaling a > Xa, where A is a nonvan- 
ishing constant. This rescaling corresponds to the coordinate rcparamctrization 
t 1— * A _1 £. Similarly to the electromagnetic case, the functional Iqr can be 
generalized to 



lGR,p[ a ] = \ 



^(a-ao) 2 + h ac h 



btl 



d t k a b\ ( d t k 



Vhd 3 x 



where here, [i > and ao is a given, strictly positive function on spacetime. The 
penalty term ^t 2 (a — ao) 2 might be useful in order to guarantee the uniqueness 
of minima and to ensure that the lapse is positive. Below, we analyze the well 
posedness of the variational principle associated to Igr,h m the weak field limit 
and show that the resulting condition for the lapse leads to an L 2 estimate for 
k a b and its first order spatial derivatives (see Lemma Q below). 

A similar functional was recently used in j.'ilj in order to construct a new 
geometric invariant measuring the amount of radiation contained in a data set 
and was also considered in |32| . Functionals very similar to Iqr also arise in 
the theory of thin elastic plates (see, for instance, |33| ) . 

2 This problem could be avoided by requiring an inhomogencous boundary condition for 
the lapse; however, for our key estimate in Lemma it turns out that we need to vary Icr 
over the space of all a £ H 2 (S). 
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2.6 The IBVP for Einstein's vacuum field equations 

With the above considerations we may formulate the following IBVP for the 
Einstein vacuum field equations. The generalization to matter fields depends 
on the precise model for the matter. For this reason, from now on, we set 

Tab = 0. 

1. Specify initial data (hij, kij) on a three-dimensional Riemannian compact 
manifold with smooth boundary dE satisfying the Hamiltonian and mo- 
mentum constraint equations (|15I16[) . 

2. Specify boundary data q(t) on [0, T] x SS. The initial and boundary data 
must satisfy certain compatibility conditions on the intersection {0} x <9E 
of the initial surface with the timelike boundary (see, for example HUE])- 

3. Choose the shift as an a priori specified vector field (3 l on [0, T] x E such 
that j3 l Ni < everywhere on T. Notice that this inequality is independent 
of the (yet unknown) metric. 

4. Find a curve t i— > (hij, kij)(t), < t < T, solution of the evolution equa- 
tions H13I14|1 . where in Eq. (|14|l the lapse is given by minimizing the 
functional Igr.u) such that 

(hij, kij) (0) = (hij, hj ) (39) 

and such that the boundary conditions 

N l M,(t) = 0, 9o(t) = c* 4 (t) + q(t), (40) 

hold, where Mi, and \&4 are defined in terms of (hij, kij) according to 
Eq. (TH and Eqs. (|32I33I35I36I) . respectively. 

In the following, we analyze the well posedness of this problem under simplifying 
assumptions. First, we only consider the weak field regime in which all the 
equations are linearized about the flat spacetime solution (hij, kij) — (5ij,0), 
a = 1, /3* = 0. Second, we assume that the domain E = (0, 1) x T 2 is a cube 
where the two opposite faces y — const and z — const are identified with each 
other (although most of the results below are valid for more general domains). 
The main result of this article is the well posedness of the above stated IBVP 
under these simplifying assumptions, see Theorem below. Results for more 
general cases will be discussed elsewhere |2*3^ . 

3 Gauge condition for the lapse 

In this section we analyze the well posedness of the variational principle for 
IcR.p,- As indicated above, for simplicity, we only consider here the case of 
linearization about flat spacetime. In this case, it is natural to set ao = 1, and 
we expand a = 1 + / and ignore all terms in Igr,^ which are cubic or higher 



12 



order in / and the perturbed metric and extrinsic curvature. This leads to the 
functional 1^ defined in Eq. I|41[) below where now Fij is equal to the linearized 
Ricci tensor. 

Before we proceed, we fix some notation: We denote by £ 2 (E, R), L 2 (E, V) 
and L 2 (E,S) the spaces of square-integrable functions, vector fields and sym- 
metric tensor fields, respectively, on E. Similarly, i? m (E, R), i? m (E, V), ff m (E, S) 
denote the Sobolev spaces of order m of smooth functions, vector fields and 
symmetric tensor fields on E and C°°(E,R), C°°(E, V), C°°(E, S) the space of 
functions, vector fields and symmetric tensor fields on E which are infinitely 
differentiable on E and such that all derivatives have a continuous extension on 
E. 

Define the bounded linear operator 

W : ff 2 (E,R) -> L 2 (E,S), / » Wtjf = DiDjf. 
We consider for each /i > and F £ £ 2 (E, S) the functional 
J M : H 2 (E,R)^R 

j ^ ^ X ^ 2/2 + ( wijf ~ f ^ {Wijf ~ Fij) ^ ^ d3x - (4i) 

Theorem 1 (Minimum of 1^) The functional J M : i? 2 (E,R) — * R has the 
following properties: 

(i) There exists a global minimum f £ i/ 2 (E,R) of I ^ . For fi > this 
minimum is unique; for (J, = i/ie minimum is unique up to the addition 
of an element in the two-dimensional set 

Af = {/ = a + a\x : a , a\ £ R}. 
(ii) A global minimum f £ i? 2 (E,M) of 1^ satisfies 

I„[f\ < \\\nh { n,sy 

(Hi) If fi = 0, a global minimum f £ i? 2 (E,R) satisfies (F— W/, Wf)L 2 Cs,s) — 
0, i.e. F — W/ and W/ are orthogonal. 

(iv) If F £ C°°(Yi,S), all global minima f £ i? 2 (E,R) automatically lie in 
C°°(E,S') and satisfy the elliptic boundary value problem 

ff + W ij (W ij f-F ij )=Q onE, 
N'NUWijf-Fij) = ondE, 
TV 1 (D 3 + V j ) (Wy f-Fij)=Q on <9E, (42) 

where T> denotes the covariant derivative on 9E induced by D. 
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Proof. The proof is based on standard arguments from elliptic theory. In a 
first step, we set X = i? 2 (£,R) and rewrite the functional 1^ in the form 



^[/] = QM(/ ! /)-2J(/) + ^||F||i 2(EiS) , 



for / € X, where the bounded bilinear form : X x X — > R and the linear 
functional J : X — > R arc defined by 

Q„(/> 9) = ^ [M 2 (/, </)^<e,r) + (W/, W.g) L2(s , s) ] , 
J(<?) = i(F,W 5 ) L2(SiS) , 
for f,g £ X. Since for all f,g £ X and all £ > we have 

M^fgLiU/j = 2(Qm(/! 5) _ JCg)) + eQ ^, g) _ (43) 

and since Qfj,(g,g) > 0, it follows that / € 1 is a global minimum of 7^ if and 
only if 

QM>9) = J(9) for all g eX. (44) 

(i) Therefore, we have to show that there exists / £ X such that Eq. I|44|l 
holds. This is done using the theorem in Appendix 1X1 which summarizes stan- 
dard results from Fredholm theory. In order to apply the theorem, we first 
notice that 

Qn(f> /) + ^ll^lltfU^.R) — oII/II-H^e,) 

for all / £ X. Next, set Z = _ff 1 (E,R) and denote by i : X — * Z the inclusion 
which is compact since S is bounded. Furthermore, we introduce the following 
linear bounded operators 



M:X^X\ f h-> Qn(f, .), 

2 



L:X^X*, f^QM>-) + kf,-)z 



R-.Z^X*, /--i(/,.) z . 

With this notation, f £ X satisfies Eq. QZty if and only if Mf — J, and 
M = L + K is the sum of the coercive operator L and the compact operator 
K = Rol. According to Theorem^! a solution to this equation exists if and only 
if J £ (kerM*)- 1 , where (M*/)(#) = (Mg)(f) for all /,p £ X and (kerM 4 )- 1 - 
denotes the annihilator of kerM*. In our case, since is symmetric, 

(Af'/Xs) = Q^gJ) = Q^f,g) = (Mf)(g) 

for all f,g £ X\ hence M* = M. If > we have kerM = {0} and so 
(kerM*) = X* . Therefore, if \i > 0, there exists a unique global minimum 
f £ X . On the other hand, if fi = 0, 

kerM = {/ e X : W tj f = 0} = {a + a r x : a , ai £ R}. 
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Furthermore, by definition, J 6 (kerM)- 1 . Hence, the existence of global min- 
ima follows from Theorem 0] 

(ii) Setting e = 1, g = — / in Eq. (|43|l we obtain 

h\f\ = ip[0] ~ Q P (/, /) < /„[0] - i||F|| 2 L2(s , s) . 
(hi) This follows from Eq. I|44|l by setting g = f '. 

(iv) This can be seen by explicitly computing the solution and using inte- 
gration by parts. 

□ 

Theorem^has the following physical interpretation: Consider the linearized 
extrinsic curvature k^ and perform an infinitesimal coordinate transformation 
Sx a I— > Sx a +X a parametrized by the vector field (X a ) = (X f , X 1 ). With respect 
to this, the linearized extrinsic curvature transforms according to 

h-> kij -Wijf, (45) 

where / = X 1 . Therefore, given k G L (E, the minimization principle 
described by 1^ provides a unique way of gauge-fixing k[ GF ^ = kij — Wij f by 
choosing / such that it minimizes 



X. 



According to Theorem ^ii), this defines for each /i>0a bounded linear map 
Pl GF) : L 2 (£, S) -> L 2 (E, S),kw k<- GF ^ whose operator nor m is less or equal 
than one. In particular, we have a unique decomposition 

hj = k\f F) +Wijf (46) 

for each k S £ 2 (E, S), where k( GF ) = P^^lc. For /i = this decomposition is 
orthogonal and Pq GF ^ is an orthogonal projector, i.e. P$ GF ^ ° P$ GF ^ = -F(j GF ^ 
and (P (Gi?) )* = Pq GF) - If k e C°°(E,5), the same holds for ^ GF ~> and it 
satisfies 

M 2 / _ W lJ k (GF) =0 on £, (47) 
N l N j k\ GF) = on as, (48) 
AT^D* + V j )k\ GF) = on 52. (49) 



One of the key properties of the gauge-fixing operator that will be 

crucial in our well posedness proof below is the following: Defining the curl and 
momentum operators by 

curl : H 1 ^^) ^ L 2 (Z,S) 

hj h-> (curlk)y ee £ kl ( i D k k l j) , 
M : H l {T,,S) -» L 2 (S, V) 

kij^(Mk)i = Dik-D j kij , 
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the gauge- fixed linearized extrinsic curvature k( GF ) = P^ GF ^k has the property 
that its H 1 norm is already controlled by its L 2 norm and the L 2 norm 3 of 
curlk( GF ) and M)<l {GF) : 

Lemma 1 Let fi > 0. There exists a constant C = > such that 

l|k|| ff i ( E,s) < C [l|k||iJ( E ,5f) + ||curlk|| La(SiS) + \\Mk\\ L 2 (sy) ] (50) 
for alike P^C 100 (£,£). 

Proof. Set (N l ) = (1,0,0) and denote by A : B, C, D, E, F indices running 
over y and z. Define the quantities 



V 



V. 



Q = 

R = 

(+) - 

A — 
(-) 

A — 

QAB = 



JVW(curlk) y - = e AM D A k Bx , 
N l (Mk) t = D x k B B -D B k xB . 



lAj 



7Aj 



--(Mky +N. l e 1]k (curlk) kl N l 
-\(MU? -NiE ijk (curlk) kl N l 



— D x k xA D A k xx , 



= D B k AB - D A k B B , 



-jAiN j s ijk (curlk) k a l B = [D x k AB - D A k xB ] 



if 



where [uab]*-^ = V AB — \5 AB 5 vcd denotes the trace-free part of a two-tensor 
v ab- Using integration by parts, we first find 



Q 2 + R 2 + 7 AB (ui +) + D A k){V B +) + D B k) d 3 x 

= J [{D x k B B f + (D A k B B )(D A k c c ) + {D x k xB ){D x k xB ) 
s 

+(D A k xB )(D A k xB )] d 3 x~2 J [(k B B )(D A k xA )]l =Q dydz. (51) 

On the other hand, using the fact that D l D^kij = fi 2 f, we find 
2 J [{D l k){Mk) t - n 2 kf] d 3 x 

E 

= 2 J {D 1 k){D l k)d 3 x~2 J [kD j k xj \l =0 dydz. (52) 

£ T 2 

Adding Eqs. lj5T|) and lf52^) together, and using the boundary conditions k xx — 0, 



3 Notice that curlk( GF ) = curlk and Mk( GF ) = Mk. 
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D j k xj + D A k xA = on 92 we obtain 



J [( J D^ s B )(Afc C c) + {p i le BB ){p i k xB ) + 2{D i k){D i k)] 



+ R l + 7 j4S (vl +) + Atfc)(V^ + D B k) + 2(D't)(Mk), - 2fkf 



r(+) 



< 



R 2 + (l + Xi) 7 AB T^ + Vi +) + K 2 (Mky{Mk) t 



d 6 x 



1 | 1 



J(n l k)(D t k)d 3 x + K 3 \\k\\l 2{ ^ s) , 



where K\ , K 2 and K3 are positive constants and where we have used Theorem 
n^ii) in the last step in order to estimate fif. By choosing K\ and K 2 large 
enough such that K^ 1 + K^ 1 < 2 it follows that we can bound the H 1 norm of 
the components k xx , k x A, k A a in terms of the right-hand side of (| 5 C)|> . 

It remains to bound the H 1 norm of k a B — ^ab — \^ab^ C c- F° r this, we 
first notice the identity 

£A C kcB = -k A c£ C B 
which is valid for any symmetric, trace- less two-tensor kAB- This implies 

2D [A k B ]c = s A b£ EF D E k FC = -£abD e k ED e D c ■ 
Using this, we find 



(D i fc BC )(D<fc BC )d 3 a: 
(D x k BC )(D x k BC ) + (D A k BC )(D B k AC ) + 2(D A k BC )(D lA k B ]c) 



d A x 



(D*k BC )(D x k BC ) + 2(D B k BU )(D /l k AC ) 



BC\ 



d A x 



1 AC 1 BD [q AB + (D A k xB ) tf ] [q C D + (D c k xD ) tf ] d 3 x 



s 

1 1' 



AC 



2V. 



(-) 



D A k l 



2V [ C ] + D c k D D 



d 3 x. 



where we have used integration by parts in the second step. This proves the 
statement of the Lemma. j— j 
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4 Well posed IBVP for the linearized vacuum 
equations 



In this section we present a well posed initial-boundary value formulation for 
the linearized Einstein vacuum equations when linearized about a flat spacetime. 
For simplicity, we also assume that S = (0,l)xT 2 although most of the results 
below hold for the general case where S C R 3 is an open bounded domain with 
C°° boundary 9E. In the following, let 

{&) = (1,0,0), K) = -J=(0,l,i), 

such that , mP , fh 3 form a complex triad adapted to the boundary. We 
assume that the linearized shift is exactly zero and fix the gauge for the linearized 
lapse by requiring that p( GF )k = k where p( GF ) = Pq GF ^ is the gauge-fixing 
projection operator introduced in the previous section. The linearized equations 
have the form 

— u = Au, (53) 
at 

where u = (h, k) denotes the components of the linearized three-metric and 
extrinsic curvature and 

A (k) = (p(^) 2 Lh)' ( 54 ) 

where Ric is the linearized Ricci operator. The idea is to represent A on an 
appropriate Hilbert space H such that A : Dh C H — ► H generates a strongly 
continuous semigroup 4 P(t) = cxp(tA) on H . The unique solution of the ab- 
stract Cauchy problem 

d 

—u = Au, t > 0, 
dt 

u(t = 0) = u , u e D H (55) 

is then given by u{t) — P(t)uo, t > 0. Furthermore, the semigroup properties 
imply the existence of constants a > 1, b € K such that ||u(t)|| < aexp(M)||uo| 
for alH > and all uq E H. In particular, the problem is well posed. 

Therefore, our task is to find a Hilbert space H and a domain Dh C H such 
that A : Dh C H — > H generates a strongly continuous semigroup. There are 
well-known sufficient and necessary conditions for this to be the case, see for 
instance |34j and |35j and references therein for generalizations to quasilinear 
operators. In our case, the Hilbert space H is motivated from our general 
considerations in Sect. and we will first define A on the space of smooth 
tensor fields (h, k) satisfying the boundary conditions (the space Dq below) . 

4 That is, P is a strongly continuous map from [0, oo) to C(H) satisfying P(0) = id// and 
the semigroup property P(t + s) = P(t) o P(s) for all t, s > 0. 
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Then, it will be shown that A is quasi-dissipative 5 and that it can be extended 
(by taking its closure in H) in such a way that the extension generates a strongly 
continuous semigroup on H. The following theorem gives the general structure 
of our Hilbert space H. 

Theorem 2 (Linear constrained evolution systems with curvature map) 

Let X, Y, Z be real Hilbert spaces, and let 

A : D(A) cl^l (main evolution vector field), (56) 
A : D{A) cF^y (constraint evolution vector field), (57) 
B : D(B) C Z — > Z (curvature evolution vector field), (58) 

be densely- defined linear operators on X , Y and Z, respectively. Further, let 
Dq C D(A) be a dense linear subspace of X and 

So : D C X — > Y ( constraint map ), (59) 
To : D C X — > Z ( curvature map ), (60) 

be closable linear operators such that Dq C D, A(Dq) C D, Sq(Dq) C D{A), 
To(Dq) C D(B) and such that the following intertwining relations hold 

S a Au = As u, (61) 
TqAu = BTqu + LiSqu + Lqu, (62) 

for all u G D where L\ : Y — > Z , L : X — > Z are bounded linear maps. 
Furthermore, assume the existence of constants do, a\, 02, a, b G M such that 

(u,Au) x < a \\u\\ 2 x +ai||S , u||y + a 2 ||T u||| for all u e D , 

(c, Ac) Y < a\\c\\\ for all c e D(A), (63) 

(w, Bw) z < b\\w\\ 2 z for all w G D{B). 

Then, 

(a) The operator Ro'.DcX^YxZ defined by R n u = (Sou,T n u) for all 
u G D, is closable. The domain D(R) of its closure R = Rq is contained in 
D(S) n D(T) where S = S~^, T = % and Ru = (Su, Tu) for all u e D{R). 
The linear space H = D(R), equipped with the scalar product 

{u,v) H = (u,v) x + {Ru,Rv) Y xz 

= (u,v) x + (Su,Sv) Y + (Tu,Tv) z , 

for u,v G H defines a Hilbert space. Furthermore, the restriction of the 
operators S and T to H are bounded linear operators from HtoY and Z, 
respectively. 

5 A is called dissipative if (u, Au) < for all u G D(A) and quasi-dissipative if there exists 
a constant b £ R such that A — b is dissipative. 
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(b) If Dq C H is dense in H, the operator Ah ■ D(Ah) C H — ► H defined by 

D(A H ) = D , A H u = Au, u e D , 
is densely- defined and quasi- dissipative. In particular, Ah is closable. 

(c) In addition, assume that (A — Ah)Dq is dense in H for some sufficiently 
large A > 0. Then, Ah generates a strongly continuous semigroup P : 
[0,oo) C(H) on H. 

(d) If in addition A generates a strongly continuous semigroup Q : [0, 00) — > 
C(Y) on Y, then 

SP(t) = Q(t) S\ H 

for all t > 0. In particular, ker S\ H is left- invariant under P(t) for all 
t>0. 

Proof, (a) Let u n be a sequence in D which converges to zero in X such that 
RoUn = {SoU n ,Tou n ) converges to some (c, w) € Y x Z . Since So and To are 
closable, it follows that c = and w = 0. Therefore, Rq is closable. Next, let 
u G D(R). Then, by definition of the closure, there exists a sequence u n in D 
which converges to u in X such that Rou n — > Ru. Since Rou n — (Sou, ToU n ), it 
follows that u G D(S) n D(T) and that Ru = (Su, Tu). Next, it is clear that H 
is linear and that (., .)h defines a scalar product on H. The completeness of H 
is an immediate consequence of the closedness of R. Finally, the fact that S\ H 
and T\ H are bounded is clear. 

(b) According to the assumptions, A(Do) C H, so Ah is a well-defined, 
linear operator. Next, let u G Dq. Using the intertwining relations I|tjll62f> and 
the estimates (|63|) we obtain 

(u, A H u) H = (u, Au) x + (Sou, S Au) Y + (T u, T Au) z 

= (u, Au) x + (Sou, AS u)y + (T u, BT u + LxSqu + L u) z 

< a \\u\\ x + (01 + a)\\S u\\ Y + (a 2 + 6)||T u||| 

+ ||Li||||Tou||z||5ott||y + ||Lo||||?bu||z||w||x 

|2 
\H 



■ KM? 



for a sufficiently large constant K > which is independent of u. Since Ah is 
densely-defined and quasi-dissipative it is closable 6 . 

(c) This is a direct consequence of the Lumer-Phillips theorem 7 . 

(d) For this, we first notice that A is closable since according to our assump- 
tions it is densely-defined and quasi-dissipative. In a next step we show that 

SA H u = ASu (64) 

for all u G D(Ah)- Let u G D(Ah)- By definition of the closure, there exists 
a sequence u n in D(Ah) = Do such that u n — > u in H and Anu n — > Ahu in 

6 See Theorem 4.5(c) in chapter 1 of Ref. 1341 for a proof. 
7 See Theorem 5.7 in Ref. 1351 for a proof. 
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H. Since for all n e N, SAnU n = ASu n and since S\ H : H — > F is bounded 

so that Su„ — > Su in y and SAnu n — > SAhu in F, it follows that 5u G P(-A) 

and that ASu = SAjju. 

Next, let u € P(A#) and define the curve c : [0, oo) — > F by c(i) = SP(t)u, 
t > 0. Since : P — > F is linear and bounded, c is differentiable on (0, oo) 
and 

4c(t) = P(t)u = I,SP(i)u = 3c(t) 
at 

for all t > 0. Since c(0) = 5u it follows by uniqueness of the Cauchy problem 
associated to A that c(t) = Q(t)Su. Therefore, SP(t)u = Q(t)Su for all u g 
£)(j4h). Since is dense in H and 5| ff : P — > F is continuous the 

statement of the theorem follows. rj 

In order to apply this theorem to the linearized vacuum equations, we define 
the operators A, A, B, So and T and the function spaces X, Y and Z as follows. 
We start with the definition of the constraint evolution vector field A. 

4.1 The constraint propagation system 

Define F = fl^E.R) x {M g P^E, V) : iV'M, = 0}, and define the operator 
A : D(A) C F -> F by 

P(A) = {(P,M) g C°°(E,R) x C°°(E, F) : JV^P = O.APMj = 0}, 
i[H\_f divM \ 
V M ) ~ V gradP J ' 

for all (P, M) g D(A), where here and in the following we use the notation 
divM = D l Mi and (gradP), = D;P. 

Lemma 2 (Well posedness of the constraint propagation system) The 

operator A is a densely-defined, linear operator on Y with the following proper- 
ties: 

(i) A is dissipative, that is, (c,Ac)y < for all c g D{A). 
(ii) For all A > 

(A-A)(P(A)) = C°°(E,R) x {M g C°°(E, V) : N l M, = 0}. 

fra) A is closable and its closure generates a strongly continuous semigroup 
Q : [0,oo) — > £(F). 

Proof. First, it is clear that A is linear and that D(A) is dense in F. 
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(i) Let c = (H, M) 6 -D(^)- Using integration by parts we find 

M * " Jl»™ tH I>H ) (D i D>M l) + M'n l H + (D>M' m D { H ) ]<i> x 
e 

= y L> 4 [ifMi + {DiH^DjMi)} (fx 

E 

= /" [^Af, + (I» 3 'flO(D J M x )]^ =0 di/tte. 

Since M x = -D x l? = the boundary integral vanishes and it follows that A is 
dissipative. 

(ii) This follows using general theorems [3] about symmetric linear operators 
with maximal dissipative boundary conditions. However, in this particular case, 
it is not difficult to give a direct proof: Let A > 0, and let F € C°°(S,R) and 
G e C°°(S, V) with WGi = 0. We want to find (H, M) e D(A) such that 

\H-D l M t = F, (65) 
AM, - DjfT = Gi . (66) 

Eqs. (|65I66|I and (H, M) S D(A) imply the Neumann boundary-value problem 

(A 2 - DW^H = XF + D i G i , 
N l D,H = 0, 

which has a unique solution H e C°°(S,R). Setting M = A -1 (G + gradi?) € 
C°°CE,V) it follows_that (i?,M) € D(A) satisfies (A - i) (H, M) = (F,G). 

(iii) Since C°°(t,R) x {M e C°°(S,T/) : = 0} is dense in Y this 
follows from (i) and (ii) and the Lumer-Phillips theorem 8 . r- j 



4.2 The Weyl propagation system 

Next, we analyze the propagation of the Weyl curvature. As discussed in Sect. 
[5] we describe the evolution of the Weyl tensor by the system constructed in pQ. 
For this, define Z = L 2 (Y,,S) x L 2 (£,S), and let B : D(B) C Z -> Z be the 
densely-defined linear operator 

D{B) = {(E,B) G C°°(t,S) x C°°(t,S) : V \E,B] = c* 4 [E,B]}, 

R ( E A _ / -(<nirlB)tf +JV (i e i) «JVj fe (divB), \ 
\B ) { (curlEjy-JV^WjVkCdivE), ^ ' 

for all (E, B) e F>(B). Here, c is a complex constant of magnitude smaller or 
equal than one and the maps (E, B) h-> W [E,B] and (E,B) h-> ^[E, B] are 

8 See Theorem 5.7 in Ref. 1351 for a proof. 
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defined by 

*o[E,B] = [E ij -e kH N k B l j ]m i m j , (67) 
* 4 [E,B] = [E ij + e k uN k B l -]m i m j . (68) 

As shown in pQ , the operator B is symmetric with maximally dissipative bound- 
ary conditions. The following lemma is a consequence of general theorems [3]. 

Lemma 3 (Well posedness of the Weyl propagation system) The oper- 
ator B is a densely-defined, linear operator on Z with the following properties: 

(i) B is dissipative. 

(li) For all A > 

(A - B)(D(B)) = C°°(E, S) x C*°°(fj, S). 

(Hi) B is closable and its closure generates a strongly continuous semigroup on 
Z. 

4.3 The main evolution system 

Next, we turn our attention to the main evolution system. Using the gauge- 
fixing projection operator p( GF ^ = Pq GF ^ introduced in the previous section we 
define the Hilbert space X = #*(£, S) x P<- GF ^L 2 (Z, S) and the dense subspace 

D = {(h, k) £ C°°(fj, S) x p( GF )C°°(£, S) : N l (Mk) l = 0}. 

In terms of the linearized Ricci operator Ric : H 2 (T,, S) — > T 2 (FJ, S) defined by 

(Rich) l3 = D {l D k h j)k - X - (D k D k h t] + 5 kl D l D,h kl ) , h e ff 2 (£, S), 

the main evolution vector field is given by A : D(A) cl^l where D{A) = D 
and 

' h \ _ / -2k 

k J - I p(GF)Rich 



A 



for all (h, k) £ D. 

Next, we define the constraint map So : D C X — > Y and the curvature map 
T : D C X -> Z by 



k J \ Mk J ' V k / V curlk / V k 

(69) 

where the linearized Hamiltonian operator H : H 2 (T,,S) — > T 2 (E,R) is defined 

by 

Hh= ls ij (Rich)ij = i (/.>7>T,, - S ij D k D k h ij ) , h e ff 2 (S,S). 



23 



Finally, we define the subspace Dq of D as the subspace of X consisting of 
smooth fields satisfying the boundary conditions. More precisely, 

D = {(h, k) e C°°(E, S) x F( GF )C°°(E, S) : 

^(Mkjj = 0, NWiiHh) = 0, * [T (h,k)] = c* 4 [T (h, k)]}, 

where V&o and "J 4 are the maps defined in Eqs. (|67I68|I and c is the complex 
constant appearing in the definition of D(B). For the proof below, we require 
that -1 which excludes the case of a "conducting boundary" described by 
m % m?Eij = 0. 

In the next Lemma and Propositions and [5] below we show that the as- 
sumptions of Theorem[2]are satisfied. This leads to the main result in Theorem 

m 

Lemma 4 (i) The operators Sq and Tq are closable. 
(11) A(D ) C D, S a (D a ) C D(A), T a (D a ) C D(B) and 

SqAu = ASqu, 

TqAu = BTqu + LiSqu, 

for all u £ Dq where L\ : Y — > Z is the bounded linear operator defined by 

T ( H \ - ( \D {i M j) -\& ii D k M k + \N^ i) ^N k ^cxa\1A)i 
V M y ^ iN^NkDtH 

for all (H,M) € F. 

(m,) There are constants ao, 0,1,0*2 € K smc/i i/iai 

(u, < Oo||«|lx + «i||5'om||^ + a 2 ||r u||| 

/or a/Z it G Do ■ 

Proof. See Appendix 151 |— I 

Following Theorem |21 we denote by R the closure of the operator Rq : D C 
X — ► y x Z defined by Rou = (Sou, T$u) for all u 6 £> and define the Hilbert 
space = D(R) equipped with the scalar product 

(u, v) H = (u, v) x + (Su, Sv) Y + (Tu, Tv) z 

for u,v € H and define the linear operator Ah ■ Do C H — > i? by 14 = Am, 
for u S £>o. It remains to prove the following two propositions: 

Proposition 1 D C H is dense. 

Proposition 2 (A — A)(D ) = D for all A > 0. 
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For the proof of Proposition 2] we rely on the following two lemmas which 
are proven in Appendix iBl 

Lemma 5 Given e > and G G C°°(dT,,R) there exists h G C°°(E, S) such 
that 

N'D^Hh) = G 

and 

\\Hw(x,s) < £, \\RicHL2(n,s) < l|Wh|| H i (s>R) < e. 

Lemma 6 Given e > and q G C°°(9S, C) i/iere existe k G C°°(E,S') smc/i 

N l (Mk), = 0, ■: curl k ;, ,//,'/„ ' = q 

and 

\\Hl*(h,s) < £) ||curlk|| L 2 (SiS ) < e, || _Mk|| H i (s>V r) < e. 

Proof of Proposition ^ Let u = (h, k) G H and e > be given. We have 
to show that there exists u G -Do such that ||u — u\\h < £• In order to do so, we 
first notice that D is dense in H and pick u = (h, k) G D such that 

e 2 

Next, we use Lemma and choose h G C°°(E, 5) such that N l Di(Hh) = — 
NWiiHh) and 

2 2 2 

l|h||ffi (s , s) < I2 ' ll ffic h|l!2 (s , S ) < Y2 ' HWh|||n (E)R) < ^2 • 

Next, we notice that the boundary condition \&o[E, B] = c^^E, B] is equivalent 
to 

■ „ 1 1 - c ■ ■ 

BjiTrfm? = Earn m? . 

3 il+c 3 

Therefore, we use Lemma[f)]and pick k G C°°(S, S) such that Aik = 0, 

1 1 - c 

curl [k + K\ijm m 3 = — [Ric (h + h)]ijm l m? 

and 

2 22 

WHh^.s) < > H curllt lli2 (SiS) < £ — , \\M\tf H ^ tV) < £ — . 

Finally, set u = (h + h,k + P< GF )k). Noticing that MP {GF) k = Mk and 
curlP^ Gp, ^k = curlk it follows that u G -Do and that 

= l|h|| 2 ffl(s , s) + ll^ (Gi;,) k||i 2(SiS) + ||Wh||^ I(E)R) + \\Mkf Hl{ ^ v) 
+ || [i«c h] TF || 1 2(SiS) + ||curlk||J a(EiS) + ttHir 
< e 2 , 
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where we have used the fact that p( GF ) is a projector in the last step. 



□ 



Proof of Proposition [5J Let A > and F G D. In the following, we 
construct u = (h, k) G Dq such that (A — A)u = P. For this, we first use 
LemmaEIii) and find (H,M) G D{A) such that 



(A -A) 



P 
M 



S Q F. 



Next, using Lemmal^ii) we find (E, B) G D(B) such that 



(X-B) 



E 
B 



T P + L x 



H 
M 



(70) 



(71) 



In a next step, we define the auxiliary variables 



Q = divE igradP, P = div B + ^curl M, 

3 2 

and prove that they vanish. Multiplying Q and P by A and using Eqs. H7()I71|) 
we obtain the system 



-\e A B D B P x 



( Ox \ 


I 


Qa 




Fx 




V Pa J 


V 



e AB D A Q B 



where A, B run over y and z. Multiplying both sides from the left with (Q x ,2Q A , P x , 2P A ) 
and integrating over S, it follows that P = Q = 0, as claimed. 
Next, we set 



F 



-2k 

p(GF) ( E + 2 5H ) 



AT 1 div E, 



Since NWiH = and Q = it follows that N' l M(E + |<5P) 2 : 
^NW.H = 0. Therefore, it follows that u G D. Next, using Eq. (JTOJ) and 
Q = we find 



A 



H 
M 



S u 



divM divTWk 





which proves that (H, M) = 5 w. Finally, Eq. P = Q = and M = Mk 

yield 



A 



E 
B 



T u 



curl B + curl k 





which shows that (E, B) = F u. Therefore, u G Do and (A — A)u = F. 
To summarize, we have shown: 



□ 
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Theorem 3 (Main result) Let c E C be such that \c\ < 1 and c ^ -1. The 

operator A : Dq C H — ► H (cf. Eq. {5$) describing the linearized Einstein 
evolution equations is closable and its closure generates a strongly continuous 
semigroup P : [0, oo) — > C(H) on H . In particular, given initial data uq = 
(h,k) G C°°(t,,S) x C°°(t,,S) satisfying the gauge condition P^^k = k and 
the compatibility conditions 

N*Di (D k D l h kl - S kl DW 3 h kl ) = 0, 
*o[T (h,k)] =c# 4 pb(h,k)], 

where 'J'g and ^4 are the maps defined in Eqs. J6'7I6'#|) and T is the curvature 
map defined in Eq. f6!A) . the curve u : [0,oo) -> ff,t h u(t) = P(t)uo is 
continuous and differ entiable on (0,oo) and is the unique solution of the abstract 
Cauchy problem l,5.5j) . Furthermore, if uq satisfies the constraints, i,e, if SqUq = 
0, the constraints are satisfied for all t > 0, i.e. Su(t) — for all t > 0. 

Remarks: 

1. For each t > the solution u(t) lies not only in H but in the domain 
D(Aji) of the closure Ajj of A^. We do not intend to give an explicit 
representation of this space in this article. 

2. It should be possible to show that the solution is smooth provided extra 
compatibility conditions are required for the initial data. 

3. Using Duhamel's principle 9 , it is not difficult to generalize the theorem to 
inhomogeneous boundary conditions and to the presence of source terms 
in the evolution equation. 

5 Conclusions 

In this article, we have discussed some new ideas for tackling the problem of 
obtaining a well posed IB VP for metric formulations of Einstein's field equations. 
These ideas go beyond casting the evolution equations into symmetric hyperbolic 
form with maximally dissipative boundary conditions which, except in some 
rather restricted situations [5J EJ |H1 El > do not seem to be flexible enough to 
be made compatible with the constraints for metric-based formulations. In 
particular, we have analyzed a gauge condition for the lapse which is obtained 
by minimizing a functional representing the norm squared of the time derivative 
of the extrinsic curvature. This leads to a fourth order elliptic boundary value 
problem. We have shown that coupling this elliptic problem to the Einstein 
evolution equations in the 3 + 1 decomposition introduced in [22] with a zero 
shift and imposing constraint-preserving boundary conditions controlling the 
Weyl scalar ^>q leads to a well posed IBVP in the linearized regime. Despite 

9 See, for example, Corollary 2.11 in chapter 4 of Ref. 1341 . 
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the fact that the results so far have only been obtained in the weak field regime, 
there are several interesting points on which we would like to comment. 

First, it is known that when the lapse is frozen or densitized, the evolution 
equations we have analyzed are only weakly hyperbolic |2I3 I3Z] ■ For such sys- 
tems there are examples of solutions with frequency-dependent exponentially 
growing modes |^] |23 E3 ■ It is therefore a priori not clear that a well posed 
IBVP can be obtained for the Einstein evolution equations in [221 ■ O n the other 
hand, as discussed in Sect. |3 even though these equations are weakly hyper- 
bolic, they induce a "nice" evolution for the constraint and curvature fields. 
The well posed initial-boundary value formulation presented in this article is 
based on this particular property of the evolution equations and on the imple- 
mentation of a fourth order elliptic gauge condition for the lapse instead of a 
frozen or densitized lapse. A well posed elliptic-hyperbolic formulation for the 
full nonlinear Einstein equations without boundaries was given in |H8j . 

The next point is related to the discretization of the problem. First, our 
formulation requires solving a fourth order elliptic equation at each (or each 
few) timesteps which might be computationally expensive. In addition to that, 
since our well posedness proof relies on the propagation of the constraints, it 
is a priori not clear that a "naive" discretization of the problem will lead to 
a stable and convergent scheme. It might be the case that one has to choose 
very special discretizations techniques such that the estimates in the continuum 
case can be mimicked at the discrete level. This may require some constraint 
projection mechanism. 

Next, we would like to remark that our approach is quite general and should 
work for any metric formulation of the Einstein equations for which the con- 
straint propagation system can be cast into symmetric hyperbolic form and for 
which sufficient regularity for the constraint fields can be shown. In particular, 
imposing the same gauge condition as in this article, we expect that a well posed 
IBVP can also be derived for families of linearized symmetric hyperbolic first or- 
der formulations [23 EH! or linearized mixed first-order second-order hyperbolic 
formulations [371 which might be more suitable for numerical discretization. 

Finally, the question remains as to what other gauge conditions may work in 
our approach. In particular, the geometrical meaning of the boundary surface 
obtained (as embedded in the spacetime constructed) must be clarified. This 
issue is likely related to the choice of the shift vector. In [2] a functional J(a, /3*) 
of the lapse a and the shift vector /3* was introduced in order to construct 
approximate Killing fields for a given data set. For zero shift our functional 
Igr[o^ is closely related to J. It should be interesting to investigate gauge 
conditions that are obtained by varying J with respect to lapse and shift. 

These questions, as well as the generalization to linearizations about more 
general spacetimes including inner excision boundaries, will be considered in 
future work. 
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A Compact perturbations of coercive operators 

In this appendix we state the following theorem which is a summary of well- 
known results. 

Theorem 4 Let X be a reflexive Banach space, and let L : X — > X* be a linear, 
bounded coercive operator, i.e. L G C(X,X*) and there exists 6 > such that 
L(u)(u) > <5||it|| 2 for all u £ X. Furthermore, let K G C(X,X*) be a compact 
linear operator and set M := L + K . Finally, let M : X — > X* denote the 
bounded linear operator defined by (M u)(v) = (Mv)(u) for all u,v G X . Then, 

(i) L : X — ► X* is invertible with bounded inverse, and < <5 _1 . 

(ii) ker M and kerM* are finite dimensional and have equal dimensions. 

(Hi) RanM — (kerM*)- 1 , where (kerM*)- 1 denotes the annihilator of ker M*. 

Proof. Let I : X — > X** denote the map defined by (Iu)(u>) — lo{u) for all 
u G X, u> G X*. Since X is reflexive, / is an isometric isomorphism. Next, 
let M* G C(X** , X*) denote the (Banach space) adjoint of M, defined by 
(M*u**)(v) = u**(Mv) for all v G X, u** G X**. Then, M t = M* o I. For the 
following, we use the formula 10 

R^A= (ker A*) 1 - (72) 

which holds for any bounded linear operator A : X — > X* with A f — A* o I on 
a reflexive Banach space X. 

(i) The coercivity of L implies that 

\\Lu\\ > 8\\u\\ (73) 

and 

||L*u|| > S\\u\\ (74) 

for all u G X. The first inequality implies that L has trivial kernel and closed 
range, the second inequality implies that L* has trivial kernel. Eq. (|72|l then 
implies that L is bijective. < S^ 1 now follows from the inequality (|73|l . 

10 See, for instance, gT] Sect. III. 
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(ii) It follows from (i) that M = L + K is Fredholm with the same index 
as L which is zero 11 . In particular, M has closed range, and dimkerM = 
dimker M l < oo. 

(iii) Since M has closed range, 

RanAf = RanM = (ker M 1 )^ . (75) 

□ 



B Proof of Lemma |4lj5] and El 

Proof of Lemma QJ (i) Suppose u n is a sequence in D which converges to 
zero in X and such that Su n converges to v £ Y. We have to show that v = 0. 
In order to see this, write u n = (h( n \ k( n )) and v = (H, M). Then, h( n ) -> in 
i? 1 (£,S , ) ) WhW -> i? in H\Y,,R) and k<") -> in L 2 (£,S), Mk^] -> M in 

Next, take a test function G Cq°(S, R). Then, 



{ip,H) L 2 (1 



lim feHhWjypj 



lim l(i?V, * W A^ ~ Dih^) LH x, v) 



= 0. 



Since this holds for all ip e C^°(S,R) it follows that H = 0. Similarly, let 

^eC °°(£,n 

W.MJi.^y) = lim (V>,A4k(")) L2(SiV) 

= o, 

which proves that M = 0. Therefore, So is closable. The proof that T is 
closable is similar. 

(ii) This follows by direct verification. 

(iii) Let u = (h,k) e D . Then, 



(u,Au) x = (h,-2k) ffl(E , s) + (k,p( GF >ffich) i2(E , s) 

1... .., 1 



< 



h|||i (s , s) + ||k||^i (SiS) + 2ll k !l^(s,s) + 2ll ffich H" 



< IN + Ci [||k||i a(EiS) + ||curlk||2 2(SiS) 

+ c 2 [\\[mch} TF \\l 2 ^ s) + \\nh\\l 2(j:j 

< a \\u\\x + ai||S u||y + a 2 ||r ii|||, 



L 2 (£,S) 
2 

L 2 (£,V) 



|Mk|' 2 



11 See, for instance, 1421. Corollaries A. 7. 2 and A. 7. 5. 
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where Ci, C2, clq, oil, 02 are constants independent of u and where we have used 
Lemma ^ in the third step. q 

Proof of Lemma [5J First, we can assume that G vanishes on one of the 
boundary components, say on 1 = 1. Otherwise, one can construct h as a 
superposition of two fields each one satisfying the statement of the lemma for 
G vanishing on x — and x = 1, respectively. In order to proceed, let "0 G 
C°°([0,oo),R) be such that 

(i) V'"(0) = 1, 

(ii) tp(x)=0 for all x > 2 

and set ip n (x) := n~ 3 ip(nx) for all n > 2 and x E [0,1]. By construction, 
tfj"'(0) = 1 and tfj n (x) = for all x S [2/n, 1]. Next, define for each n > 2 

4° (a, y, 2) = (2^iV, - 5 lJ )Mx)G(y, z). 
Since ^"(«) = ip"'(nx) = for all x E [2/n, 1] it follows that 

l|hW|| 2 ff3(s>s) <^ 

for a sufficiently large constant if independent of n. Therefore, the £f 3 -norm of 
h^™) can be made arbitrarily small. Furthermore, one finds 



WhW = -ip"{nx)G{y, 



so NWiCHh^'^l = G. 



□ 



Proof of Lemma |5J Similarly to the proof of Lcmma[5]we may assume that 
q vanishes oni= 1. Next, let \& E C°°([0, 00), R) be such that 

(i) tf(0) = 0, 

(ii) *'(0) = 1 

(iii) $(x) = for all x > 2 

and set x P„(x) := n^ 1 'i>(nx) for all n > 2 and x E [0,1]. By construction, 
<jr„(0) = 0, *;,(0) = 1 and * n (as) = for all a; G [2/n, 1]. Next, define for each 
n > 2 



fcjy (a;, y, z) = -i^ n (x)q(y, z)m k mi 



c.c. 



where c.c. denotes the complex conjugate of the previous expression. Since 
*&' n (x) = ^>'(nx) = for all x E [2/n, 1] it follows immediately that there is a 
constant K > such that 

||kW||?, 1(S)S) <f 
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for all n > 2. Therefore, the i/^-norm of k^™) can be made arbitrarily small. 
Next, since ty n depends on x only one finds 

(Mk {n) )i = i^ n fhifh k D k q + c.c, 

hence Mk^\ x=0 = and ||.Mk( n )|| 

l 2 (s,v) can be made arbitrarily small. Fi- 
nally, 

(curlk (n) ) fc; = ^\nx)m k rhiq + ^ n {x)N( k rh l) rh :l Djq + c.c. 
which implies that (curl k^)/~im k m l | = q. rj 
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